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Abstract 



Fuchsian equations provide a way of constructing large classes of space- 
times whose singularities can be described in detail. In some of the ap- 
plications of this technique only the analytic case could be handled up 
^ ■ to now. This paper develops a method of removing the undesirable hy- 

' pothesis of analyticity. This is applied to the specific case of the Gowdy 

, spacetimes in order to show that analogues of the results known in the an- 

' alytic case hold in the smooth case. As far as possible the likely strengths 

\l ' and weaknesses of the method as applied to more general problems are 

. displayed. 

o 

O '. 1 Introduction 

^ , The theory of Fuchsian equations has been applied to analyse singularities in 

W)' a variety of classes of spacetimes in general relativity. The term 'Fuchsian 

equations' has not always been used in the literature on this subject and in this 
paper it denotes a certain class of singular differential equations in a generic way. 
The existing results in this area will be surveyed briefly below. In this approach 
' spacetimes containing singularities are parametrized by some functions which 

play the role of data on the singularity. In some cases it was necessary to assume 
the analyticity of these functions. In other cases smoothness was sufficient. The 
aim of this paper is to develop ways of removing the analyticity requirement. 
These will be illustrated by the case of Gowdy spacetimes which represent an 
ideal laboratory for testing new ideas on the mathematical treatment of the 
Einstein equations. 

It may not be immediately obvious why the apparently technical distinction 
between analytic (C"^) functions and smooth (C°°) functions should be signifi- 
cant with a view to physical applications. There are at least two reasons why 
it is important. The first is that the physical notion of causality cannot be 
reasonably formulated within the class of analytic functions, since the unique 
continuation property of the latter means that the solution of an equation at 
one point influences its behaviour at all other points, and not only at causally 
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related points. Connected with this is the fact that solutions of an equation do 
not depend continuously on initial data in any useful sense. For more discussion 
of these points see 0, particularly section 2.4. The second, which is also a direct 
consequence of the unique continuation property of analytic functions, is that 
there is not the same freedom to construct solutions with certain interesting 
properties within the analytic class. An example of this will be given in section 
^ below. 

There are a number of results on Fuchsian equations with smooth coefficients 
in the literature and to start with we need to understand why these do not 
apply directly to Gowdy spacetimes. The general form of a system of Fuchsian 
equations for a vector-valued unknown function u is 

Qu 

t—+N{x)u^tf{t,x,u,D^u) (1) 

Here a; is a point in some Euclidean space and D^u is a shorthand for the first 
order derivatives of u with respect to the spatial variables x. The function / is 
required to extend continuously to i = while the matrix N{x) is required to 
satisfy some positivity condition, which may depend on the particular context. 
An example of a condition of this type is that the eigenvalues of the matrix 
N{x) should have non-negative real parts for all x. Note that the apparently 
more general system 

Oil 

t—+N{x)u^ef{t,x,u,D,u) (2) 

with a > can be reduced to the form ([^) by introducing t" as a new time 
variable. This results in the matrix N{x) being rescaled by a factor a, but this 
does not affect its positivity properties. 

One approach to proving existence theorems for Fuchsian systems which 
does not require any analyticity assumption is due to Claudel and Newman |^ . 
Of course, in the context of data which arc merely smooth (or even finitely 
differentiable) the system must be hyperbolic. This is needed to prove existence 
for an equation without any singular behaviour in t and the singularity cannot 
be expected to improve the situation. In the Claudel-Newman theorem it is 
required that / have an asymptotic expansion about t — in integral powers 
of t and this property is inherited by the solution. The positivity assumption 
on N{x) is weaker than that mentioned above. It is only required that there be 
no eigenvalues which arc negative integers. An important element of the proof 
of the theorem is to expand the candidate solution u in powers of t, writing 
it in the form uq + ui where uq is an appropriate polynomial in t and mi is a 
remainder of higher order in t. Then ui solves a Fuchsian system where the 
eigenvalues of N{x) are shifted by an integer in the positive direction. If it is 
possible to expand to a sufficiently high order then the shifted eigenvalues all 
have positive real parts. The condition which makes the expansion possible is 
that a polynomial uq can be found which satisfies the equation up to an error 
of sufhciently high order in t. The only obstruction to this is if the shifted 
eigenvalue becomes zero at some stage in the process and this is prevented by 
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the assumption made on the eigenvalues of the original matrix. In cases such 
as the system arising in the analysis of Gowdy spacetimes in the solutions 
cannot be expanded in integral powers of t. Instead non-integral (and even 
x-dependent) powers of t and logarithms occur. For this reason the method of 

does not apply directly. One of the main methods of the present paper is to 
extend the technique of using expansions of the solution to shift the eigenvalues 
of N{x) to cases where terms more complicated than integral powers of t occur. 

Another approach to proving existence theorems for Fuchsian systems with 
smooth coefficients is due to Kichenassamy Q . In that case there is no restriction 
that solutions have expansions in integral powers of t. It is, however, required 
that the matrix N(x) be independent of x. This is not satisfied in the Gowdy 
case. Since the eigenvalues of the matrix N{x) correspond to powers occurring 
in the expansions, this is a consequence of the dependence of these powers on x 
in the system coming from the Gowdy spacetimes. Thus the result of |^ does 
not apply. It might be possible to extend the proofs in ||^ to the case of non- 
constant N, but this will not be attempted here. Both the proofs of Claudel and 
Newman and of Kichenassamy involve the use of sophisticated techniques from 
functional analysis, namely semigroup theory and the Yosida approximation, 
respectively. These will be avoided in the approach developed in the following. 

In the paper [|j Anguige treats the case of plane symmetric solutions of the 
Einstein equations coupled to a perfect fluid. He uses a direct energy argument 
of a type familiar in the theory of regular symmetric hyperbolic systems. It is 
important in his proof that N(x) is independent of x and positive semi-definite. 
This kind of argument will be generalized in the following. On the one hand the 
condition that N{x) should be constant will be removed. On the other hand the 
positivity condition on the matrix will be relaxed by a method related to that 
used in . The eigenvalue condition used in ||] will be replaced by a condition 
of formal solvability which abstracts its essential significance. 

Next a brief survey of the literature on applications of Fuchsian equations 
to general relativity will be given. It appears that the first paper containing 
an application of this kind was where Moncrief proved the existence of 
a large class of analytic spacetimes with analytic compact Cauchy horizons. 
Later Newman |6| based his work on isotropic singularities and Penrose's 
Weyl curvature hypothesis on existence theorems for hyperbolic systems with 
singularities of Fuchsian type. These papers did not include proofs of the re- 
quired theorems but the necessary proofs were provided in More recently, 
results on isotropic singularities for more general matter models were obtained 
by Anguige and Tod [|[ ^ |^. Their theorems require no symmetry assump- 
tions but are confined to a special type of singularities. On the other hand 
Anguige 1^ proved a theorem on the existence of non-isotropic singularities in 
plane symmetric spacetimes with perfect fluid. 

Another line of development of the applications of Fuchsian equations in 
general relativity starts with the work of Kichenassamy and Rendall |jri| on sin- 
gularities in analytic Gowdy spacetimes. It builds on previous work of Kichenas- 
samy outside general relativity (see [|lOj and references therein). This direction 
is continued in the papers B, M and which deal with vacuum models 
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with two spacelike Killing vectors, solutions of the Einstein-scalar field equa- 
tions and analogues of Gowdy models in string cosmology, respectively. In all 
these cases analyticity is assumed. A notable feature of the result of is that 
it makes no symmetry assumptions and so, on the basis of function-counting 
arguments, concerns general solutions of the Einstein equations with the given 
matter model. 

The paper is organized as follows. In the second section the notion of a 
formal solution of a Fuchsian system is defined. Assumptions on the coefficients 
are exhibited which guarantee the existence of a formal solution. They are 
fulfilled by the first order form of the Gowdy equations introduced in In 
the third section this form of the equations is modified slightly so as to obtain a 
symmetric hyperbolic system. Its formal solvability is shown to follow from that 
of the original system. The system satisfied by the remainder term which is the 
difference between a true solution of the system and an approximate solution 
is computed. The fourth section proves an existence theorem which is general 
enough to apply to the case of Gowdy spacetimes with sufficiently low velocity 
{k < 3/4). In the fifth section yet another form of the equations is used to cover 
the remainder of the full low velocity case (fc < 1). The wider applicablity of 
the methods of the paper is discussed in the final section. 



2 Formal solutions 

If the function / in (||) is smooth at t = and hence admits an asymptotic 
expansion about t = in integral powers of t then it can be useful to expand 
it in this way. In the following a generalization to less smooth functions / will 
play an important role. This uses the notion of a formal solution of equation 
which will now be defined. It will also be convenient for the following to 
introduce a notion of regularity of functions adapted to the given situation. An 
analogous notion in an analytic context was introduced in 
Definition 1 A function z{t, x) defined on an open subset of [0, oo) x i?^ and 
taking values in is called regular if it is C°° for all i > and if its partial 
derivatives (defined for t > 0) of any order with respect to the variables x G 
extend continuously to t = 0. 

Definition 2 A finite sequence (wi, U2, ■ ■ ■ , Up) of functions defined on an open 
subset U of [0, oo) x R" containing {0} x R" is called a formal solution of order 
P of (0) if 

1. each Ui is regular 

2. tdtUi + N{x)ui — tf{t, X, Ui, DxUi) — 0{t^) for alH as i 

Here and in the following the 0-symbols are taken in the sense of uniform 
convergence on compact subsets. 

In jl^ an iteration was defined which, in the case that the function / has 
suitable analyticity properties, converges to a solution of (|l|). It is doubtful if 
it converges in any useful sense when / is only smooth. However it can be used 
to produce a formal solution of any desired order, as will now be shown. 
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Lemma 2.1 If the function / is regular and the matrix N{x) is smooth and 
satisfies an estimate of the form ||cr^^^^)|| < C with a constant C for ct in a 
neighbourhood of zero then for each p the equation (|l|) has a formal solution of 
order p which vanishes at t = 0. 

Proof First some definitions from will be recalled. For a regular function 
u define F[u] = tf{t, x, u, D^u). Then F[u] is also regular and is 0{t) as t — > 
together with all its spatial derivatives. If v is regular and 0{t) as f — > 
together with all its spatial derivatives, define u = H[v\ = Jp a^^^^~^v{at)d(7. 
Then u is regular and 0(t) together with all its spatial derivatives and satisfies 
{tdt + N)u — V. Then, if G is defined to be the composition FH, any fixed point 
V oi G defines a solution u of (|l|) hy u = H[v]. Let ui = 0. It defines a formal 
solution of of order one. It will be shown that defining Ui = HG^~^[ui] 
defines a formal solution of order p for each p. Note the relation Ui+i = HFui. 

The first defining property of a formal solution is easily proved by induction. 
The main point is to verify the second property. To do this it will first be shown 
that Ui+i — Ui = 0{t'^) for each i, and that spatial derivatives of all orders of 
the Ui satisfy the corresponding estimates. For i = 1 this follows directly from 
the properties already demonstrated. To prove the general case, consider the 
equation obtained by forming the difference of the equations satisfied by Ui+i 
and Ui. This gives 

tdt{ui+i - Ui) + N{x)(ui+i - Ui) = tMi(ui - Ui^i) + tM2Dx{u.i - Ui^i) (3) 

for regular functions Mi and M2 of the arguments t, a;, m^, D^Ui, DxUi-i, 

obtained by applying the mean value theorem to differences. The right hand 
side of (|^) is 0{V-). Then the fact can be applied that the operator H preserves 
the set of functions which are 0{t^ ) for any j . Spatial derivatives can be handled 
in the same way. To complete the proof of the lemma, consider the relation: 

tdtUi+i + N{x)ui+i - tf{t, X, Ui+i,DxUi+i) 

-t[f{t, x, Ui+i, D^Ui+i) - f(t, X, Ui, D^Ui)] (4) 

Using the mean value theorem and the estimates already obtained for Ui^i — Ui 
shows that the right hand side of (^) is 0(<*+^). 

Remark A general criterion for checking the condition on N(x) required to 
apply this lemma has been given in (|^). 

Next some basic equations for the Gowdy spacetimes will be recalled. More 
details can be found in The basic unknowns are two real- valued functions 
X{t, x) and Z{t, x) of two variables. New variables u and v are defined so that 
the relations 

Z(t, x) = k(x) log t + 4){x) + t^u{t, x) (5) 

and 

X{t, x) = Xoix) + {^p{x) + v{t, x)) (6) 

hold, where k, Xq, 4> and ip are given functions. The positive constant e will be 
restricted later. Next introduce further variables by setting 

[uo,Ui,U2,Vo,Vi,V2) = {u,tdtU,tUx,V,tdtV,tVx) (7) 
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The Gowdy equations imply the following first order system: 



tdtui 



Ui 

-2eui - e^uo + t^~^{Kx logi + (t>xx) + tdxU2 
- exp(-20 - 2t^uo){t^^~''{vi + 2fcwo + ^kipf 
-2t^^''Xox{v2 + tipx + kx{vo + i^)tlogt) 
-t^''-'-{v2 + tiPx + 2kx{vQ + i:)t\ogtf} 
tdx{uo + ui) 



(8) 



t 



_2-2k-e 




tdtU2 

tdtvo 



(9) 
(10) 

(11) 



tdtvi = -2kvi+t'^-'^''Xoxx+tdxiv2 + tijx) + 4:kxiv2+t^x)tlogt 
+{vo + i^)[2kxxt^ logi + 4{kxt\ogt)^] 
+2f{vi + 2kvo + 2kip){ui + euo) 
-2Xoxt'-^^{kx \ogt + 4,x+ t'dxUo) 

-2t{dx{vo 2kx{vo + V-) \0gt){kxtl0gt + t(j)x + t'w2) (12) 

tdtV2 = tdx{vo + vi) (13) 

Here some minor errors in the equations given in [Tl| have been corrected^. This 
system is of the form (^) which implies, as indicated in the introduction, that it 
can be brought into the form (|l|) by a change of time variable. This possibility 
will be used freely without further comment in the following. After the change 
of time coordinate the system arising from the Gowdy equations satisfies the 
hypotheses of Lemma 2.1 provided e < 2fc and e < 2 — 2k. In particular, the 
bound on a^'^^^ was verified directly in ||ll|. Alternatively, it follows easily from 
the criterion given in . Hence the above system has a formal solution of any 
order which vanishes at t = 0. 

3 The symmetric hyperbolic system 

In the following the form of the Gowdy equations introduced in section ^ will be 
called the first reduced system. Now it will be modified to get a system which, 
while less convenient for showing formal solvability, is symmetric hyperbolic 
and therefore appropriate for allowing the theory of hyperbolic equations to be 
applied. It is obtained from the first reduced system by making the substitutions 
U2 = tdxUo and V2 — tdxVo in some places. The result is: 

tdtUQ = ui (14) 



tdtUi = -2eui- e^UQ+t^ "(kxx'^ogt + (j)xx) + tdxU2 

-exp(-2(?!)- 2i'Mo){^^''"'(wi + 2fcwo + 2kijj)^ ~i 
-2t^-''Xox{v2+ti^x + kx{vo + ^P)tlogt) 
-f'''^V2+t^Px + 2kx{v„+^)t\ogt)^} 



.2-2fc-e 




(15) 



thank Aurore Cabet for pointing these out 



6 



tdtU2 

tdtvo 



U2 + tdxUi 



(16) 
(17) 



tdtvi = -2kvi + <2-2feXcte^ + td.^(v2 + + 4fc^(v2 + ti;^)t\ogt 
+{vo + ip)[2k^xt^ \ogt + A{kxt\ogtf] 
+2f^{vi + 2kvo + 2kip){ui + euo) 
-2Xo^t^-^^{k^ \ogt + (j)^ + t'-^U2) 

-2{v2 + tdxi> + 2tk:^{vQ + ip) logi)(A;^Hogt + tcj)^ + i'u2) (18) 
tdtV2 = V2+ tdxVi (19) 

This system, which will be referred to as the second reduced system, has the 
advantage of being symmetric hyperbolic but also has two potential disadvan- 
tages. The one is that the matrix N{x) has been modified while the second 
is that possibly dangerous powers of t have been introduced on the right hand 
side. The matrix N{x) acquires two negative eigenvalues, which will have to be 
dealt with by appropriate methods in due course. As far as the other problem 
is concerned, the power ^i+'^-^fc introduced. This power should be positive. 
This can be achieved subject to the inequalities already assumed for e if and 
only if A; < 3/4. This restriction appears unnatural, but will be assumed in this 
section and the next for the second reduced system. 

In fact, although the desired inequalities relating e and k(x) can be ensm'cd 
by a suitable choice of e at any given point x, this cannot be done at all points 
simultaneously by a single choice of the constant e. It can be ensured in a 
neighbourhood of any given point. Solutions will be constructed in local neigh- 
bourhoods of this kind and then put together using the domain of dependence 
to get a solution which is global in x. 

Any formal solution {ui, . . . , Up} of the first reduced system which vanishes 
at t = is also a formal solution of the second reduced system. For tdt{{u2 — 
tdxUa)i) ~ 0{t^) and using the fact that the formal solution vanishes a,t t = 
it follows that {u2 — tdxUo)i = 0{f). This can then be used to see that 
the difference terms arising when passing from the first to the second reduced 
systems are 0{t^), assuming the condition k < 3/4. Although it is not of 
significance for the following it is interesting to note that for 3/4 < < 1 the 
sequence whose element with index i is the element of the formal solution of the 
first reduced system with index i + 1 is a formal solution of the second reduced 
system. 

Given a formal solution it is possible to consider the difference between an 
actual solution and the formal solution and the equations which this difference 
satisfies. The hope is that this equation is more tractable analytically than the 
first and second reduced systems. This is a generalization of the procedure of 
subtracting a Taylor polynomial of finite order in the case that the solutions 
are smooth at t = 0. Rather than doing this calculation in the specific case 
of the Gowdy system it will be done for the following more general symmetric 
hyperbolic Fuchsian system: 



tdfU + N{x)u + tA' {t, x, u)djU = tf{t, x, u) 



(20) 
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The condition for symmetric hyberbolicity is that the matrices should be 
symmetric. As before, all coefficients in the equation are assumed regular. If 
{ui} is a formal solution and Vi — t^^{u — Ui) then 

tdtVi + {N{x) +il)vi +tA^t,x,u, +fvi)djVi = tgi{t,x,Vi) (21) 

for some regular function gi. Choosing i large enough ensures that the eigen- 
values of N + il have positive real parts, or even that the matrix is positive 
definite. In the Gowdy case this will be referred to as the third reduced system. 

4 The basic existence theorem 

In this section a local existence theorem will be proved for solutions of the 
Gowdy equations in a neighbourhood of the initial singularity in the case that 
the data k, X^, and ■0 are merely smooth. In the case of analytic data the 
problem was solved in JlTf . Thus if the data are approximated by a sequence of 
analytic data {km,XQm, 4>m, V'm), a corresponding sequence of analytic solutions 
is obtained. At the same time formal solutions can be obtained for both the 
approximate data and the actual data. Denote the former by Umi, where the first 
index corresponds to the sequence of data and the second to the enumeration of 
elements of an approximate solution. Denote the latter by u^, as before. If the 
approximate solutions are constructed as in the proof of the Lemma 2.1 then 
Umi Ui as m — > cxD, uniformly on compact subsets. The same is true of the 
spatial derivatives of these functions of any order. It can be concluded that the 
sequence of coefficients obtained for the third reduced system is also convergent 
on compact sets as m — > oo. The sequence of solutions of these equations which 
we have is only defined a priori on a time interval which depends on m. However, 
using the global existence theorem for the Gowdy equations |jl^ it is possible to 
conclude that a sequence of smooth solutions exists on a common time interval. 
The aim now is to show that this is a Cauchy sequence in a sufficiently strong 
topology. If that can be done then it will follow that the sequence converges to a 
limit which is a solution corresponding to the smooth data originally prescribed. 

The tool to obtain convergence of the approximations is the technique of 
energy estimates. This requires some preliminary remarks on linear algebra. 
Consider the matrix-valued function N{x) in the Fuchsian system. Spatial 
derivatives of a solution of this system also satisfy a system of the same form, 
but with a different matrix N{x). Suppose, for instance, we consider the first 
derivative D^u of the unknown. The system for the pair (u, D^u) has a matrix 
in its singular term which has diagonal blocks N{x) and an off-diagonal block 
involving Dj.N{x). This off-diagonal block does not affect the eigenvalues of 
the matrix but may well affect whether it is positive definite or not. Since the 
positive definiteness of matrices like this is important in what follows we adopt 
a strategy which avoids positivity being lost. In order to implement this it will 
be assumed that N{x) is positive definite. Then use the variable w = KD^u for 
a positive constant K instead of D^-u itself. For the equation satisfied by {u, w) 
the matrix of interest is positive definite provided K is chosen sufficiently small. 



8 



The same trick works for higher derivatives. It suffices to replace the collection 
of unknowns {D"m} by — K^°'^D°'u. Let the matrix corresponding to N 
arising in the system for all these derivatives up to order s be denoted by N^-^\ 
By construction it is positive definite. 

The standard method of energy estimates (see e.g. Q or, for a discussion 
aimed at relativists, |^]) proceeds by estimating the Sobolev norms of solu- 
tions. The usual Sobolev norm is given by HwUff^ = {'^\a\<s ll^^^llia)^^^- For 
the present purposes it is convenient to use the equivalent norm ||u||//s.a' = 
(X]|a|<s II^^^IIl^)^''^: where K is, as before, a small enough positive con- 
stant. 

As a first application of energy estimates, a theorem on the domain of de- 
pendence will be proved. Let u and v be two regular solutions of ( po|) vanishing 
at t = 0. Then their difference satisfies an equation of the following form: 

u — v) + A' {t, X, u)dj{u — v) = M {t, x){u — v) (22) 

where M is a regular function constructed from u and v. Choose two times ti 
and t2 with < ti < t2 and let G be the region defined by the inequalities 
ti < t < t2 and I a; I < 2^2 — t. Let S'l and ^2 be its intersections with t = 
ti and t = t2 respectively. Now multiply the equation (^ ) by e~^'^{u — v), 
integrate over G and integrate by parts in the way this is usually done in the 
derivation of energy estimates. Here k is a positive constant. The singular term 
containing N can be discarded, due to its sign, giving an estimate of the form 
g-nt2j^ < e~''*i/i -I- Ig, where Ii and I2 are the norms of the restrictions 
of M — w to 5i and ^2 respectively, and Ig is a volume integral over G which 
for K sufficiently large is negative unless u — v is identically zero on G. Letting 
ti tend to zero, so that Ii —>■ 0, gives a contradiction unless u — v ~ on G. 
Thus it has been proved that the solutions u and v agree on G. This proves a 
domain of dependence property for solutions of (^0|) which can be used for the 
purpose of gluing together solutions. This means that even if we are interested 
in producing solutions on manifolds it is enough to solve the problem on R". 
Moreover, it is possible to consider without loss of generality the case of cut-off 
coefficients and data. By this we mean that there is a compact subset of R" 
such that for x outside this compact set the initial data vanish, the coefficients 
and / vanish and N is constant. 

The aim is now to construct solutions of the third reduced system for Gowdy 
in the case of smooth coefficients. As already indicated it is enough to do this 
under the assumption of a cut-off in space. A sequence of functions which is a 
candidate for a sequence converging to a solution of the third reduced system 
has already been produced. This sequence is obtained by fixing a value of i 
sufficiently large that the matrix occurring in the singular term of the third 
reduced system is positive definite and forming the difference of the solution of 
the second reduced system corresponding to the data {km, 4>in, ^Om, V'm) and the 
function Umi- It will now be shown by using energy estimates that this sequence 
is bounded in suitable Sobolev spaces and in fact is a Cauchy sequence. To do 
this differentiate with respect to t and substitute the third reduced 
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system into the result. This gives 

d/dt{\\u\\H^^Kf) = -t-'{N^^^u^^\u^^^)L2^j, + R (23) 

Here u^'*^ is the cohection of all derivatives of u up order s and R is the sum of the 
terms which arise in the regular case, i.e. in the case where N is identically zero. 
The first term on the right hand side is non-positive and may be discarded. The 
terms in R can be estimated just as in the regular case and this gives a bound 
for the iJ** norm of u, provided s > n/2 + 1. Next the difference of successive 
approximants will be estimated. An attempt to apply the standard techniques 
in order to show that the sequence is Cauchy only meets one difficulty not 
present in the regular Cauchy problem. This is due to differences of the matrices 
N for successive elements of the sequence. This can be overcome in a way 
similar to that used above where the norms were scaled. The trick is to consider 
the collection of the derivatives of the functions of the sequence up to order s 
together with the derivatives of the differences of successive elements of the 
sequence up to order s — 1, the derivatives of the differences being multiplied by 
an additional factor K. This once again ensures that the final matrix obtained 
is positive definite. 

It follows from the above discussions that the sequence of solutions Vi of the 
third reduced system converges to a solution of the third reduced system corre- 
sponding to the original smooth data. This can then be used to define solutions 
of the second and first reduced systems and finally a solution of the Gowdy sys- 
tem itself corresponding to the data we started with. The construction is such 
that the interval on which convergence is obtained may depend on s. However 
standard results about symmetric hyperbolic systems with smooth coefficients 
show that the solutions for all values of s can be extended to a common time 
interval. Hence a solution is obtained on that time interval which is C°° for 
t > 0. The results of this discussion are summed up in the following theorem, 
which may be compared with Theorem 1 of 

Theorem 4.1 Let k{x), Xq{x), (j){x) and ip{x) be C°° and assume that < 
k{x) < 3/4 for all x. Then there exists a solution of the Gowdy equations with 
following properties. For each spatial point x there exists an open neighbour- 
hood Ux of X and a number > such that the restriction of the solution to 
Ux satisfies (^) and (|^) with e — Ex, where u and v are regular and tend to zero 
as t — > 0. The Ux and ex can be chosen in such a way that the inequalities 
2k{y) ~ 1 < ex < min{2fc(y),2 — 2k{y)} are satisfied for all x and all y G Ux- 
Under these conditions the solution is unique. 

Remark A formulation of the theorem which is equivalent but cleaner can be 
obtained by replacing the constant e by a function e{x). Then it would not be 
necessary to introduce the Ux ■ 

In Q a theorem was proved concerning high velocity analytic solutions in 
the case where Xq is independent of x. The method used to prove Theorem 4.1 
applies straightforwardly to the high velocity case to give an analogue of the 
result of in the smooth case. The following theorem results. 
Theorem 4.2 Let k{x), (j){x) and ip{x) be C°° and let Xq be a constant. As- 
sume that k{x) > 0. Then there exists a solution of the Gowdy equations with 
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the following properties. For each spatial point x there exists an open neigh- 
bourhood Ux of X and a number > such that the restriction of the solution 
to Ux satisfies and with t — ex, where u and v are regular and tend to 
zero as t ^ 0. The Ux and Cx can be chosen in such a way that the inequality 
ex < '2k{y) is satisfied for all x and all y G Ux- Under these conditions the 
solution is unique. 

5 Data with intermediate velocity 

In the previous section an existence theorem was proved for Gowdy spacetimes 
under the restriction 0</c<3/4on the function k. Next it will be shown that 
using a different ansatz allows the range 1/2 < k < Itobe treated. The two 
together then cover the whole range < k < 1 for which results were available 
in the analytic case. The new ansatz involves expanding the function Z to a 
higher order in t. The ansatz for X remains unchanged. Now Z is of the form 

Z^klogt + (j) + at^-^'' + t^-^'^+'u (24) 

where a = (2 — 2k)~'^X^^. Reexpressing the Gowdy equations in terms of the 
new variables u and v and reducing to first order as before leads to an analogue 
of the second reduced system of section ^. The evolution equations for uq, U2, vq 
and V2 are the same as before and will not be repeated. The modified equation 
for Ml is 

tdtui = -2{2-2k + e)ui~{2~2k + efuo + t'^''^''{kxx^ogt + ct)xx)+tdxU2 
+t^-''[taxx -4:kx logt{tax + t'u2) + 4:klt{\ogtfia + t'uo)] 
- exp(-2(/) - 2at2-2fc „ 2t^-^''+'ua){t'^''-^-''{vi + 2kvo + 2kipf 
-^2t^''-^-'Xox{v2 + tijjx + kx{vo + i^)tlogt) 

-t'^^'^-'{v2 + ti^x + 2kx{vQ + i!)t\ogtf} (25) 
and that for vi is 

tdtvi = -2kvi + t^-^'^Xoxx + tdx{v2 + txjjx) + Akx{v2 + tTjjx)t\0gt 
+ {vo + ip)[2kxxt^ logt + 4{kxtlogt)^] 
+2t^-^''+'{vi + 2kvo + 2ki}}){ui + (2 - 2fc + e)uo) 
+ (4-4fc)i2-2fe(i,^ +2fct;o + 2fcV)a (26) 
-2Xoxt'-''\kx \ogt + (j)x + t^-^'^ax + t^-^''+'u2 
-2fc^t^"2fc \ogt{a + t^UQ)) ~2{v2 + tdxijj + 2tkxivo +tp)\ogt) x 
{kxt\ogt + t(f>x+t^^^''ax+t'^^^''^'u2 - 2kxt^^^''logt{a + fuo) 

There is also an obvious analogue of the first reduced system of section ^. The 
existence of formal solutions of the latter is guaranteed by Lemma 2.1 and these 
give rise to formal solutions of the second reduced system as in section ^ 

Consider now the sequence of analytic solutions of the Gowdy equations 
corresponding to a sequence of analytic approximations to the smooth data 
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of interest. It will be shown that, under the condition that 1/2 < fc < 1, 
these define a sequence of regular solutions of the second reduced system of this 
section. To do this it is necessary to show that for each of these solutions the 
function Z admits an asymptotic expansion of the form ( p^ ) and not just of the 
form (|^), which is known a priori. To do this it suffices to apply the existence 
theorem of Q to the first reduced system of this section with the analytic data. 

Once these facts are known, it is straightforward to prove an existence the- 
orem for the second reduced system of this section using the same techniques 
as were applied to the second reduced system of section || provided certain in- 
equalities are satisfied. These are the inequalities which ensure that each term 
on the right hand side of the equation contains a positive power of t. Under the 
assumption that 1/2 < A; < 1 this can be achieved by choosing the positive real 
number e to satisfy 4fc — 3<e<2fc — 1. The following theorem is obtained. 
Theorem 5.1 Let Xo(x), and tli{x) be C°° and assume that 1/2 < 

k(x) < 1 for all x. Then there exists a solution of the Gowdy equations with 
the following properties. For each spatial point x there exists an open neigh- 
bourhood Ux of X and a number e^; > such that the restriction of the solution 
to Ux satisfies (|^ and (||) with e = e^, where u and v are regular and tend to 
zero as t 0. The Ux and ex can be chosen in such a way that the inequalities 
4:k{y) — 3 < Ex < 2k{y) — 1 are satisfied for all x and all y G Ux- Under these 
conditions the solution is unique. 



6 Discussion 

The theorems stated in this paper have all concerned Gowdy spacetimes. It is 
nevertheless clear that many of the arguments are much more generally applica- 
ble. At the same time some steps are essentially related to the specific Gowdy 
case. A general discussion of the procedure will now be given which separates 
the general from the particular as much as possible. The first step is to make a 
suitable ansatz for the solutions to be constructed as the sum of an explicit part 
and a remainder. There may be more than one useful way of doing this. For 
example, in the Gowdy case equations (1) and (|) were useful for proving one 



theorem while replacing (|5|) by ( |24[ ) allowed a different theorem to be proved. 
The second step is to reduce the equations to first order. The aim is to produce 
a system of Fuchsian form for which the theorem of |^ ensures the existence of 
solutions corresponding to the case where the free functions in the ansatz are 
analytic. If these free functions are merely smooth the lemma proved in section 
2 may be used to show the existence of formal solutions. 

The third step is to produce a system which is symmetric hyperbolic and in 
Fuchsian form. At this stage the matrix N{x) may have negative eigenvalues, 
as is the case in the Gowdy example. It needs to be shown that the formal 
solutions already produced define formal solutions of the symmetric hyperbolic 
system. From this point on the argument proceeds on a general level, with no 
more details of the Gowdy special case being used. 

It is instructive at this stage to consider what difficulties would be likely 
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to arise in an attempt to generalize the results of from the analytic to the 
smooth case. One problem is to bring the equations into a suitable hyperbolic 
form by the choice of coordinate or gauge conditions. There was no difficulty of 
this type in the Gowdy case, where a rather rigid preferred coordinate system 
is available. In more general cases it will be necessary to choose a form of the 
reduced Einstein equations carefully from the myriad on offer. If a symmetric 
hyperbolic system is obtained it is likely to involve a matrix A" multiplying the 
time derivative of the unknown which is not the identity, thus going beyond the 
case discussed above. Even worse, it may be difficult to ensure that remains 
bounded and uniformly positive definite as i ^ 0. These conditions are very 
important for the use of energy estimates. 

To conclude the paper, an application will be presented where the flexibil- 
ity of smooth functions is essential. Existence theorems have been proved for 
Gowdy spacetimes in the low velocity case and, under the condition that Xq is 
constant, also in the high velocity case. Using the domain of dependence these 
can be combined to give a more general class of solutions. To do this consider 
a smooth function Xq which is constant on a non-empty open interval /. Now 
complete this to data {k,Xo, 0, tp) in such a way that fc < 1 on the closure of the 
complement of /. Then each point x has a neighbourhood on which one of the 
existence theorems applies and the resulting local solutions can be put together 
to produce a solution corresponding to the chosen initial data globally in x. If 
we tried to do this construction with analytic data then Xq would have to be 
globally constant and nothing new would be obtained. 
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